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TECHNICAL PUBLICATION 


AN ESTIMATION OF THE LIKELIHOOD OF SIGNIFICANT ERUPTIONS 
DURING 2000-2009 USING POISSON STATISTICS ON TWO-POINT 
MOVING AVERAGES OF THE VOLCANIC TIME SERIES 


1. INTRODUCTION 


Volcanoes are found on every continent of the Earth and undersea. According to Sinikin and Siebert, 1 
more than 1,500 volcanoes have been active in the past 10,000 yr, accounting for -8,000 known eruptions. 
About 400 volcanoes eaipted during the 20th century. Some two-thirds of the volcanoes are located in the 
Northern Hemisphere and, similarly, some two-thirds are found to lie around the Pacific Ocean margin, 
reflecting a predilection of the occurrence of volcanic eruptions to tectonic plate movements. 

Episodes of temporary cooling have long been suspected of being related to the occurrences of 
large volcanic eruptions. For example, Franklin 2 suggested that perhaps the eaiption of an Icelandic volcano 

was responsible for the unusual weather of 1783 1784 in Europe, and Lamb 3 described many of the unusual 

meteorological effects that followed particular large volcanic eruptions. These effects, lasting for months 
to years, include strange sunsets and reduction of the direct measurement of normal-incidence solar radiation 
by 10 percent or more, in addition to temporary cooling. 4-11 Furthermore, large volcanic eruptions have 
been linked to enhanced sulfate concentrations in ice cores 12-14 and to changes in tree ring densities. 15-1 7 

According to Robock, 11 large volcanic eruptions influence climate by injecting large amounts of 
chemically and microphysically active gases and solid aerosol particles directly into the stratosphere. 18-21 
These gases and aerosol particles have an e -folding residence time of ~1 yr. They affect the Earth’s radiation 
balance and disturb the chemical equilibrium of the stratosphere, leading to cooling at the Earth’s surface 
(0. 1-0.2 °C) and heating in the stratosphere, with these effects often lasting 2-3 yr or more. Following the 
eruption, the cloud of sulfate gases and aerosol particles spread outward from the eruption site, encircling 
the globe and expanding toward both poles in a matter of weeks. The actual degree of spreading depends 
upon the particular distribution of winds at the time of the eruption and the latitude of the volcano. The 
sulfur species (mainly sulfur dioxide) react with the hydroxy group (OH) and water to form sulfuric acid. 
It is the resulting sulfuric acid aerosols that produce the dominant radiative effects. 

Recently, Hyde and Crowley 22 noted that a number of climate model predictions indicate that the 
average global temperature should rise about 0. 1 -0.3 °C over the next 1 0-1 5 yr. Furthermore, they suggested 
that the occurrences of large volcanic eruptions might mask the global warming effects of carbon dioxide. 
Specifically, on the basis of a 600-yr chronology of Northern Hemispheric eruptions inferred from ice 
cores, they determined probability density functions that allowed for the estimate of the probable occurrence 
within the next 1 0 yr of one or two climatically significant volcanic eruptions. They found that the probability 
of a radiative perturbation of -1 W nr 2 or larger occurring in the next 10 yr is about 35-40 percent, the 



probability of two such eruptions is ~1 5 percent, and the probability of three such eruptions is ~5 percent. 
(They found the probability of no such eruptions to be ~44 percent.) Also, they deduced that the probability 
of a Pinatubo-scale eruption (volcanic explosivity index (VET) = 6) in the next decade is ~20 percent, and 
the probability for an El Chichon-scale eruption (VEI = 5) in the next decade is ~25 percent. Hyde and 
Crowley conclude “that there is a reasonable chance of a volcanic eruption influencing the evolution of the 
anthropogenic signal over the next 10 years.” 

In this study, the number of known cataclysmic volcanic eruptions per decade over the past 
250 yr (between 1750—1 993), based on the compilation of Simkin and Siebert, 1 is determined. A cataclysmic 
volcanic eruption is one having a VEI >4, and it is these that are believed most likely to influence climate. 
(Strictly speaking, it is the sulfur richness of the eruption that determines whether or not it will have a 
significant effect on global climate.) A simple two-point moving average is applied upon the resultant time 
series to ascertain the expected number of cataclysmic volcanic eruptions during the present decade (2000- 
2009). Then by means of Poisson statistics, various probabilities of occurrence are estimated. Thus this 
study provides an alternate way of determining the likelihood that a climatically significant volcanic eruption 
will occur within the next 10 yr, as compared to that of Hyde and Crowley. 22 
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2. RESULTS 


Figure 1 plots the number of cataclysmic volcanic eruptions per decade, known to have occurred 
between 1750 and 1993, taken directly from Simkin and Siebert. 1 The eruptions are divided by hemispheric 
subgroupings: Panel (a) is the extra-tropical Southern Hemisphere (N(ETSH)), panel (b) is the tropics 
(N(T)), panel (c) is the extra-tropical Northern Hemisphere (N(ETNH)), and panel (d) is the combined data 
set (N(Combined)) which ignores the hemispheric location of the eruption. The number of events, mean, 
and standard deviation are given for each subgrouping. Additionally, for the combined grouping, a two- 
point moving average is displayed (the heavy line), showing that prior to about 1860, the number of 
cataclysmic eruptions per decade appears fewer than thereafter — an indication, perhaps, that the record 
might be less reliable before the mid- 1 800’s . Across the top are filled triangles, denoting the occurrences of 
volcanic events with a VEI >5. These particular events are identified for convenience in table 1. Such 
events are routinely identified as being among the largest of stratospheric-aerosol-producing volcanic 
eruptions. 7 ’ 11 ’ 13 ’ 15 

Using the combined data set, one finds that since 1750, there have been at least 123 cataclysmic 
volcanic eruptions. Slightly more than half occurred in the tropics and only -4 percent occurred in the 
extra-tropical Southern Hemisphere. The largest value yet observed for the two-point moving averages is 
found to have occurred in the 1910’s, having a value of 8.00 (being associated with the peak number of 
actual events for the entire sample which measures 1 1 .00). Since then, the two-point, moving averages have 
remained relatively fiat, varying between 5.00 and 6.00, at least, until the 1980’s when it suddenly jumped 
to 6.50, (The trend of the two-point moving averages since the 1930’s appears to be directed upward.) 


Figure 2 displays the variation of the first differences (Id’s) of the two-point moving averages. The 
first difference is computed as the next two-point moving average minus the present two-point moving 
average. For 16 of 23 decades (70 percent), the observed fd has been within one unit of the last available 
value of the two-point moving average. Hence, one approximates the next two-point moving average as 
being equal to the last available two-point moving average plus/minus one unit. Thus the decade of the 
1990’s should have a two-point moving average of -6.50 + 1.00, providing it is not a statistical outlier, 
implying that the actual rate of cataclysmic volcanic eruptions for the present decade (2000 -2009) should 
be -7 ± 4. 

Table 2 gives the observed decadal rate (r) of cataclysmic volcanic eruptions and the associated 
Poisson probability distribution for each of the groupings shown in figure 1 . The Poisson distribution is 
quite useful for measuring the probability of occurrence in positively (rightward) skewed discrete 
distributions, where knowledge of the mean number of random events that occur within a given time 
interval is plainly known. 23 One finds that at least two events per decade have been seen since 1750, using 
the combined data set. The probability of having no cataclysmic volcanic eruptions during the decade of 
2000-2009 is computed to be only 0.7 percent, implying that the probability of having at least one event is 
>99 percent. 
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The number of occurrences per decade of cataclysmic volcanic eruptions 
(1750’s-2000’s). The decade marked 1800’s includes one event that has no 
VEI value associated with it, but instead was identified on the basis of an ice 
core reading and tree ring density pattern — the 1809 “unknown” eruption. 
(Confirmed by Robock, private communication.) See text for details. 


Table 1. Volcanic eruptions having VEI >5 (1750-1993). 


Year 

Volcano 

Latitude 

VEI 

1800 

St. Helens 

46.2°N 

5 

1815 

Tambora 

8.3°S 

7 

1822 

Galunggung 

7.3°S 

5 

1835 

Cosiguina 

13.0°N 

5 

1853 

Chikurachki 

50.3°l\l 

5? 

1854 

Sheveluch 

56.7°N 

5 

1883 

Krakatau 

6.1 °S 

6 

1886 

Okataina (Tarawera) 

38.1 °S 

5 

1902 

Santa Maria 

1 4.8°N 

6? 

1907 

Ksudach 

51 ,8°N 

5 

1912 

Novarupta 

58.3°N 

6 

1932 

Azul, Cerro (Quzapu) 

35.7°S 

5+ 

1956 

Bezymianny 

56.0°N 

5 

1980 

St. Helens 

46.2°N 

5 

1982 

El Chichon 

17.4°N 

5 

1991 

Pinatubo 

15.1°N 

6 

1991 

Hudson, Cerro 

45.9°N 

5 


'Taken from Simkin and Siebert 1 



Decades 

Figure 2. The first difference of the two-point moving average, based on the combined 
data set. See text for details. 


From table 2, one finds that the actual observed range of cataclysmic eruptions per decade is 2-1 1 . 
Presuming that the decadal rate (r) will be 7 ± 4 for the decade of 2000-2009, from the Poisson probability 
distribution one computes the probability for such a rate to be -87 percent. Hence, one infers a very strong 
likelihood that several cataclysmic eruptions will be seen during this decade. 

Table 3 gives the observed decadal rate (r) of VEI events >5 and the associated Poisson probability 
distribution for the combined data set. The probability of occurrence of one of these largest of events within 
the next decade, with the associated large stratospheric-aerosol injections, is ~49 percent. The probability 
of at least two such events is ~15 percent, while the probability is ~4 percent that at least three such events 
will occur. Although not shown, the probability of at least one VEI event >6 occurring within the next 10 yr 
is =1 8 percent and the probability of at least two such events occurring is ~2 percent. It is only =0.1 percent 
for three or more such events occurring. 
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Table 2. The distribution of the decadal rate (r) of cataclysmic volcanic 
eruptions (1750-1993) and the associated Poisson probability 
distribution for selected groupings of eruptions. 




umumui 

HBI 

r 

EHH 

EBI 

mu 


mu 

| 

mm 

P{r) 

0 

20 

0.8187 

4 

0.0805 

3 

0.1108 

0 

0.0073 

1 

5 

0.1637 

4 

0.2028 

5 

0.2438 

0 

0.0359 

2 

0 

0.0164 

6 

0.2555 

6 

0.2681 

5 

0.0883 

3 

0 

0.0011 

5 

0.2146 

6 

0.1966 

3 

0.1449 

4 

0 

0.0001 

2 

0.1352 

5 

0.1082 

3 

0.1782 

5 

0 

0.0000 

2 

0.0681 

0 

0.0476 

4 

0.1754 

6 

0 

0.0000 

1 

0.0286 

0 

0.0174 

5 

0.1438 

7 

0 

0.0000 

0 

0.0103 

0 

0.0055 

2 


8 

0 

0.0000 

1 

0.0032 

0 

0.0015 

1 


9 

0 

0.0000 

0 

0.0009 

0 

0.0004 

1 


10 

0 

0.0000 

0 

0.0002 

0 

0.0001 

0 


11 

0 

0.0000 

0 

0.0001 

0 

0.0000 

1 


m 

0.20 

2.52 

2.20 

4.92 


Table 3. The distribution of the decadal rale (r) of VEI >5 eruptions 
and the associated Poisson probability distribution. 


r 

N(Combined) 

Obs 

P{r) 

0 

13 

0.5066 

1 

7 

0.3445 

2 

5 

0.1171 

3 

0 

0.0265 

4 

0 

0.0045 

5 

0 

0.0006 

6 

0 

0.0001 

m 

0.68 


6 













3. DISCUSSION AND SUMMARY 


Since the 1880’s, a prominent feature of long-term temperature studies has been the appearance of 
warming. 24 ' 30 This warming is directly attributable to both natural climate variability and the increasing 
effects of human activity on climate. 31-34 Differentiating the relative proportion of each has been and 
continues to be both scientifically and politically important. 35 

Recently, Karl et al 36 have reported that during the interval of May 1997 to September 1998, each 
month broke the previous monthly all-time-record high temperature. In fact, researchers have found that 
the inferred average temperature in the 20th century is probably the warmest of the past five centuries. 37 

As previously noted, Hyde and Crowley 22 have recently examined the statistical properties of a 
proxy record of climatically significant volcanic eruptions to determine the likelihood of future eruptions 
within the next 10—15 yr. Furthermore, they noted that the occurrence of one or more major volcanic 
eruptions likely would confuse the debate regarding temperature trends and the effects of human activity 
on climate. Their results indicated that there is indeed a reasonably good chance that a major volcanic 
eaiption will occur within the next 10 yr, one that might influence the evolution of the anthropogenic 
signal. Specifically, they determined that there is a 20-percent probability for a Pinatubo-scale eruption 
(VEI = 6) and a 25-percent chance for an El Chichon-scale eruption (VEI = 5), both archetypes known to 
have temporarily cooled the Earth’s surface air temperatures. 7 - 1 138 

This complementary study, using actual counts of cataclysmic volcanic eruptions (i.e., those of VEI 
>4) over the past 250 yr, taken from Simkin and Siebert, 1 and Poisson statistics, determined the probability 
distribution for such events. Over the past 25 decades, each has seen at least two cataclysmic volcanic 
eruptions, usually located in the tropics or extra-tropical Northern Hemisphere. On average, nearly five 
such eruptions per decade were seen, with the largest number having occurred in the 1910’s (1 1 events). 
Many of these cataclysmic events are known to have produced short-term global cooling, as evinced from 
surface air temperature records and proxies (i.e., ice cores and tree rings). The actu al proportion of cataclysmic 
events that have known associations with short-term global cooling arguably is ~0.1 to nearly 1, depending 
upon the size, location, and time of year of the eruption. For the largest of events (VEI >5), there appears to 
be an ~49 percent chance that at least one will occur within the next 10 yr and an ~15 percent chance that 
at least two will occur. The usual behavior of the first difference of the two-point moving averages of the 
number of cataclysmic volcanic eruptions suggests that one should expect ~7 ± 4 during the present decade, 
having a probability of occurrence of ~87 percent. The probability of at least one El Chichon-scale event 
(VEI == 5) Is computed to be 47 percent (slightly higher than that found by Hyde and Crowley 22 ), while It 
is ~ 14 percent that at least two will occur and =3 percent that at least three will occur. The probability of at 
least one Pinatubo-scale event (VEI = 6) is computed to be 15 percent (slightly less than that found by 
Hyde and Crowley 22 ), while there is an ~1 percent chance that at least two will occur, and <0. 1 percent that 
at least three will occur. The probability of at least one event where VEI is >6 is 18 percent. Therefore, as 
previously recognized by Hyde and Crowley 22 on the basis of proxy volcanic data (sulfate measurements 
in the ice core record), it seems quite reasonable to expect at least one climatically significant volcanic 
eruption within the next 10 yr. 
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